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We study the twist defects in the toric code model introduced by Bombin [Phys. Rev. Lett. 
105, 030403 (2010)]. Using a generalized 2D Jordan-Wigner transformation, we show explicitly 
the twist defects carry unpaired Majorana zero modes. We also draw the same conclusion using 
two alternative approaches of a perturbation theory and a projective construction. We propose a 
quantum non-demolition measurement scheme of the parity of Majorana modes. Such a scheme 
provides an alternative avenue to demonstrate the non-Abelian statistics of Majorana fermions with 
measurement-based braidings. 


I. INTRODUCTION 

Quantum systems exhibiting topological order have 
stimulated a lot of excitement in the past few decades be¬ 
cause a new paradigm beyond Landau’s symmetry break¬ 
ing theory is needed to describe those state^®^. In par¬ 
ticular, topology rather than symmetry should be used to 
characterize topological states. The existence of topolog¬ 
ical degeneracy of ground states makes topological sys¬ 
tems a promising candidate for topological quantum com¬ 
putation (TQC)PE1 as such degeneracy depends only on 
the topology and hence is robust against local perturba¬ 
tions. 

One important example of topological order is the toric 
code modeJ^ which is exactly solvable. With local in¬ 
teractions between spins, the system has topologically 
degenerate ground states with long-range entanglement. 
Local perturbations only introduce exponentially small 
splittings of the ground state degeneracy, and the ground 
state subspace is thus protected from local probes and 
possesses a long-range topological order. Anyons emerge 
as excitations in the toric code model. Because fusion of 
two such anyons has a definite outcome, these are abelian 
anyons and have no computational power if we want to 
encode quantum information in the fusion channels. 

An interesting variation of the toric code model that 
supports non-Abelian anyons was proposed by BombirP 
via introducing twist defects into the model. By verify¬ 
ing the fusion and braiding rules, the twists were shown 
to mimic the behavior of more exotic non-Abelian Ising 
anyons which can be used to realize the quantum gates of 
the Clifford group for TQCP. An important distinction 
was made by You and WeiP that those twists are not in¬ 
trinsic non-Abelian anyons themselves since they are not 
excitations of the Hamiltonian. Instead, twists are ex¬ 
trinsic defects with projective non-Abelian statistic^^. 
With a mean-field treatment twists are shown to carry 
unpaired Majorana fermion^. More recently, the topo¬ 
logical entanglement entropy was calculated to show that 
twist defects have the same quantum dimension and fu¬ 
sion rules as Ising anyonJ^. The twists have also been 
proposed as qubits for quantum computation in a surface 
code implementation to reduce the space-time costl^. 

In this work, we explicitly show the emergence of un¬ 


paired Majorana zero modes associated with the twist 
defects using three different appro aches: 1) a generalized 
2D Jordan-Wigner transfer mat iorP^^ 2) a projective 
construction from a Majorana surface code model; 3) a 
perturbation theory from the Kitaev honeycomb lattice 
modefSK^. In addition, we show that the parity opera¬ 
tor of two unpaired Majorana modes is indeed the logical 
operator defined in the spin language, which can be used 
to manipulate the degenerate states of the correspond¬ 
ing twists. To demonstrate the non-Abelian statistics of 
twists through braiding, one needs to move the twists by 
changing the HamiltoniarP. Here, instead we want to 
realize the braiding of Majorana fermions without actu¬ 
ally moving them. We investigate the measurement-only 
topological q uantu m computation scheme using forced 
measurementand find that the measurement-based 
braiding can be simulated efficiently using a single cycle 
of topological charge measurements followed by a sin¬ 
gle logical qubit Pauli operation. This way, we avoid 
the uncertainty associate with forced measurement whose 
number of measurements is probabilistically determined. 
Even through signatures of Majorana fermions have been 
identified in several spin-orbit coupled materials placed 
close to super conduct or^^M^, demonstrating the non- 
Abelian statistics of Majorana fermions in the condensed 
matter setting is still a very challenging task as it re¬ 
quires a combination of experimental capabilities includ¬ 
ing braiding, fusing and measuring two Majorana zero 
modes. We instead show that it is possible to verify 
the non-Abelian statistics of Majorana fermions using 
the twist defects in the surface code modeP^lt^ on a 
2D planar lattice. We believe all the required opera¬ 
tions are within experimental reach in the surface code 
setup given recent significant experimental progress in su¬ 
perconducting circuit^^^H^. Therefore, our scheme pro¬ 
vides an alternative approach to demonstrate the non- 
Abelian statistics of Majorana fermions. If supplemented 
with a si ngle- qubit tt/S gate through the magic state 
distillatiorP^, the measurement-based braiding of Ma¬ 
jorana fermions can be used for universal quantum com¬ 
putation in the surface code settin^^. 

Here is the outline of the paper. In Section|H| we first 
introduce the toric code model with twist defects and 
show explicitly there are unpaired Majorana fermions as¬ 
sociated with the twists. In Section [ni| we show how 
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to demonstrate the non-Abelian statistics of Majorana 
fermions using the twist defects. In particular, we give 
a detailed account of the measurement-based braiding of 
Majorana fermions and show that it can be done in a sin¬ 
gle cycle of parity measurements. Finally, we conclude 
in SectionlVl 


II. UNPAIRED MAJORANA ZERO MODES 
FROM TWIST DEFECTS 


Because of their exotic non-Abelian statistics, Majo¬ 
rana fermion^^ have generated tremendous interest in 
the condensed matter community and a huge amount 
of effort has been devoted to the search of Majorana 
fermions in the past few yearP^H^. So far, majority of 
the activity has focused on topological superconductors 
supporting localized zero-energy modes^^^^. However, it 
has been shown by Petrova et that introducing 

topological defects in the Abelian phase of the Kitaev 
honeycomb modefSl provides an alternative way to gen¬ 
erate unpaired Majorana zero modes. Previously, twist 
defects in the toric code model are shown to exhibit Ising 
anyon like behavior through either fusion and braiding 
ruie^BHini or a mean-field treat menfP. Here, it is our aim 
to demonstrate explicitly the emergence of unpaired Ma¬ 
jorana zero modes associated with the twist defect s as 
has been done in the Kitaev honeycomb modeP^^. 

The toric code model was initially introduced by 
Kitae\i^ and later reformulated by WerP. To be self- 
contained, here we describe the fundamental aspects of 
the toric code model. The model is defined on a 2D 
lattice with spin-1/2 particles seating on each site, see 
Fig.B We start with the reformulated Hamiltonian with 
twist defects introduced by BombirP: 

Htcm = - ^ ( 1 ) 

k 

where the plaquette operator Aj^ is a four-body inter¬ 
action term between spins living on the corners of the 
plaquette k and is given by 


^k = 


2 

1 


Z X 
X z 


3 

4 


= A1Z2A3Z4. 


( 2 ) 


Here, X’s and Z’s are the spin-1/2 Pauli operators. For 
the plaquette containing a twist defect, the plaquette op¬ 
erator has to be modified to 


4 = 


= XiZ2X3Z^Y5, 


(3) 


where Y = iXZ. 

It is easy to check that all the plaquette operators com¬ 
mute with each other: each Ak (or A^) shares two spins 


with its neighboring plaquettes and anti-commutation re¬ 
lations between different Pauli operators guarantee A^ 
(or A'j^) commutes with them; for all the other plaque¬ 
ttes, Ak (or A^) simply shares zero spin with them. This 
nice property makes the toric code model exactly solv¬ 
able. Since A| = (A'^)^ = 1, each plaquette operator 
has two eigenvalues ±1. Therefore, the ground state en¬ 
ergy corresponds to simultaneously minimizing the con¬ 
tribution from each plaquette. This leads to the condi¬ 
tion that each Ak (or A'^) takes the eigenvalue +1 and 
the ground state is a common eigenstate of all the pla¬ 
quette operators: Ak\GS) = |G5') for all k. When the 
lattice is placed on a torus, the ground state has a de¬ 
generacy of 4 X for an even x even lattice, and 

2 X for an even x odd or odd x odd lattic^. 

The factor of 4 or 2 is the number of distinct set of exci¬ 
tations allowed by the underlying lattice topology. 
and Apiaq denote the number of sites and plaquettes, 
respectively. In contrast, with an open boundary condi¬ 
tion, the ground states support gapless edge modes due 
to Majorana fermion^. 

Excitations can be generated by applying spin opera¬ 
tors, say Z on site 1 in Eq. It flips the states of Ak 
together with A/^’s diagonal neighbor that shares site 1 
with Ak. Hence, excitations are localized on the plaque¬ 
ttes and they are always produced in pairs with a min¬ 
imal energy gap of 2. To label the excitations, one can 
color the plaquettes into dark and light groupJ^. This 
labeling can be done consistently on the lattice for the 
regions free of twist defects. Then a electric charge e 
(magnetic charge m) is attached to excitations living in 
a dark (light) plaquette. Both e and m can be moved 
by applying spin operators to flip the plaquette states, 
for instance, we can move the charge on Ak in Eq. ([^ 
by applying Z on site 3. Connecting those spin flipping 
operators gives a string operator having its ends on two 
plaquettes where the excitations are generated. Once we 
keep moving one of the charge until the two ends of the 
string operator meet on the same plaquette again, the 
two charges will then annihilate and the system is back 
to the ground state. This tells us that fusing two e or m 
charges gives vacuum 1. The composition of one e charge 
and one m charge however gives a fermion e = exm. The 
fermionic nature of e can be verified by exchanging a com¬ 
posite pair of e X m with another pair of e x m, which 
gives rise to a 7r-phase shift in the wave-function. 1 , e, m 
and e together forms the quantum double model of zM 
They are abelian anyons because fusing any two of them 
always has a definite outcome. Eor an even x even lat¬ 
tice placed on the torus, these four types of charges can 
be distinguished and the 4-fold degeneracy of the ground 
state due to the lattice topology. While for an even x odd 
or odd X odd lattice, e becomes m if moved across the 
boundary and hence they are no longer distinguishable 
any more, so are 1 and e, giving the 2-fold topological 
degeneracy of the ground stat^. 

Introducing twist defects provides an alternative ap¬ 
proach to change the ground state degeneracy, n twist 
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FIG. 1. (color online) Planar surface code lattice with an 
open boundary condition. The spins live on the corners of 
the square lattice. The dashed blue line shows the ordering 
of the spin along the Jordan-Wigner path. The green ovals 
show one possible representation of the parity operator of the 
two twist defects. After multipling the plaquette operators 
in the shaded region, the parity operator becomes the one 
showing in red circles. 


defects will reduce the number of plaquette A^piaq by n/2 
and hence increase the degeneracy by . Each twist 
has a quantum dimension of \f2 which coincides with 
that of Major ana fermions. This connection to Major ana 
fermions has generated quite some interestPH^ mainly be¬ 
cause of the generation of non-Abelian anyons from an 
Abelian phase. Of course, caution has to be applied be¬ 
cause twist defects are not intrinsic anyons because they 
are not part of the excitation spectrum of the Hamilto¬ 
nian. Instead, they are extrinsic anyons with projective 
non-Abelian statistic^^. But still, twist defects have 
been shown to have some computational poweJ^ and can 
be used to encode qubits for quantum computatiorPl. 
Here, we show explicitly that twists are indeed associated 
with unpaired Majorana zero modes using three different 
techniques. 

In this paper, we assume a 2D planar lattice with an 
open boundary condition. Since we are only concerned 
with the topological degeneracy introduced by the twist 
defects (instead of the degeneracy associated with lattice 
topology), the boundary condition becomes irrelevant. 
More importantly, a 2D planar lattice is more realistic 
from the experimental point of view. Hereafter, we will 
refer to the 2D planar model as the (planar) surface code 
model if we are considering a Hamiltonian or surface code 
if we are using a stabilizer formulation without Hamilto¬ 
nian. 


A. Jordan-Wigner transformation of the surface 
code model 


To define a Jordan-Wigner transformation on a 2D lat¬ 
tice, we order the spins in a quasi-lD manner as shown in 
Fig-i Following the path, we define the following map¬ 


ping between the spin operators and Majorana fermions 

y = = y = (4) 

where Uj = ^ non-local disorder operator 

attached to guarantee that Pauli operators at different 
sites commute. Here, j labels the position of the spin 
along the Jordan-Wigner path. / < j simply means the 
spin at / is ordered before the spin at j on the path. 

and are Majorana fermion operators non-local in 
terms of the spin operators 

= (5) 

j'<j 

( 6 ) 

o'<j 

and they satisfy the anti-commutation relationship 

{ 7 “'> 7 “} = {if ,lj} = > ilj'^lj} = 0 - ( 7 ) 

Now, we substitute Eq. Q into the square plaquette 
operator defined in Eq. (§ and obtain 


= 


> 


Z X 
X z 




= (C7i7t)(C7272“)(C7373)(77474“)- (8) 


Depending on the direction of the Jordan-Wigner path, 
there are two possible mappings of Ak. If the path goes 
from site 1 to 4, and 3 to 2 in a counter-clockwise fashion. 


then Eq. ^ becomes 


= (n?74")(^73"72")- 

(9) 

Otherwise, we have 


4fe,0 = (77b4)(*7372)- 

(10) 


Similarly, the pentagon plaquette operator defined in 
Eq. can be rewritten as 


4 = 


Z X ' 
X Y 


which reduces to 


A' 

- 

A' 

^k,0 


= (f7i7i) (77272) (77373) (77474) (*75 75), 

( 11 ) 

= (ni“74“)(773“72“), (12) 

= (*7b4)(*7372)- (13) 


From Eqs. it is clear that each plaquette oper¬ 

ator consists of the product of two pairwise link terms 
of a-type (b-type) Majoranas when the path ordering is 
counter-clockwise (clockwise). For each site except the 
twist defects [site 5 in Eq. (pT|], both the a-type and b- 
type Majorana modes are involved in the plaquette oper¬ 
ators because the top two plaquettes containing the site 
have an opposite path ordering of the bottom two plaque¬ 
ttes. However, the twist defects are special because the 
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Majorana modes associated with them do not show up in 
the pentagon plaquette operator. Since the assignment 
of a-type ( 7 ^) and b-type ( 7 ^) Majorana fermions is ar¬ 
bitrary in Eq. 0 . the two path orderings are equivalent 
upon exchanging 7 ^ with 7 ^. Therefore, it is sufficient 
to assume the case of clockwise ordering in Eq. (13). In 
this case, 75 are involved in the two plaquettes below 
and 75 becomes an unpaired Majorana (it does not 
appear in the Hamiltonian.) Clearly, introducing a twist 
defect gives rise to an unpaired Majorana mode in the 
new representation. 

A pair of twists will produce two unpaired Majorana 
zero modes which can be used to encode qubit. In Eig.[^ 
they correspond to 7^2 and 7^9 The parity operator as¬ 
sociated with these two unpaired Majorana zero modes 
is given by 


P — '^7i27i9* (1"^) 

Now, we want to express P in the spin representation. To 
eliminate the contribution from the boundary, we need 
to modify the form of the boundary spin operators in the 
Jordan-Wigner transformation. Specifically, at sites 14 
and 15 (Eig.Q, we exchange X 14 with ¥ 14 ^ and Z 15 with 
Tis in Eqs. Q-(|^. As a result, X 14 and Z 15 replace 
Yi 4 and Tis in the definition of the disorder operators 

Ujyis = (ni<i4 l^i)V4^15(ni>fe>15n). 

Using the modified Jordan-Wigner transformation 
from Eqs. ©-(HI, the parity operator can be rewritten 
as 


P = iZi2Yi3Xi4Zi,YieYuYisXig. (15) 

We highlight the spin representation of P in green ovals in 
Fig-E Since each plaquette operator has an eigenvalue 
of -|-1 in the ground state, we can multiple P by the 
plaquette operators shaded in dark brown in Eig.[^ We 
end up with an equivalent parity operator shown in red 
circles in Eig.[^ 

P' = iXi2YiiYioZoX2oZio^ (16) 

Compared to P, P' is a better representation because it 
does not go through the boundary and has a distance of 
roughly the separation of the two twist defects. This is 
also the reason why we modify the Jordan-Wigner trans¬ 
formation; otherwise, we will have a term of Z 14 X 15 from 
the boundary in the definition of P'. It is fairly easy to 
check that P' commutes with all the plaquette operators 
and it can not be written as the product of plaquette op¬ 
erators. In fact, if we encode information in the degen¬ 
eracy states due to the two twists, then P' is essentially 
the logical Z operator accessing those states. 


B. Projective construction from a Majorana 
surface code model 

Inspired by the projective construction due to Wen^^, 
here we start with a surface code model of Majorana 



FIG. 2. (color online) A Majorana fermion model on a 2D 
lattice with twist defects. Each site has four Majorana modes 
(dots): a, 5, c, and d. a and c modes form bonds (red) in the 
horizontal direction, while b and d modes form bonds (green) 
in the vertical direction 


fermions, show there are unpaired Majorana modes at 
the twist defects, and then map it to the surface code 
model of spins via projection. The model is defined on 
a 2 D planar square lattice with four Majorana modes 
^a,b,c,d on each site n. Each site forms a pairwise 

bond of Majorana modes with its nearest neighbors in 
the form shown in Eig.|^ Basically, 7 " and 7 *^ are paired 
along horizontal bonds, and 7 ^ and 7 ^ are paired along 
vertical bonds. By construction, the Hamiltonian takes 
the form 

^MF = —^A/e, (17) 

k 

where Ak is plaquette operator of Majorana modes. Eor 
a square plaquette, it reads 



Eor a pentagon plaquette with a twist defect, it reads 





= (* 7 ^ 2 ") (nl75“) (*75^7?) 


(19) 

This model is exactly solvable because each link opera¬ 
tor = ^ 7 ^ 7 ^ or = ^ 7 ^ 7 ^ commutes with the 
Hamiltonian. Since each link operator has two eigenval¬ 
ues of ± 1 , the ground state is the common eigenstate of 
all A/c operators with eigenvalue + 1 . 

Note that in Eq. (19), 75 is absent from the pentagon 
plaquette operator and it is not involved in the two square 
plaquettes below A^ either. Therefore, 7 ^ is an unpaired 
Majorana mode associated with the twist defect. To per¬ 
form the mapping to spins, we first define two fermion 
operators out of the four Majorana modes at each site 




7n + 


— 


7n + hn 


( 20 ) 





















5 




FIG. 3. (color online) Perturbative limit of honeycomb lattice 
model with twist defects, (a) The Kitaev honeycomb model 
with two twist dislocations. The unpaired Majorana modes 
are highlighted by the dots, (b) In the strong z-bond limit 
(think green bonds), under perturbation theory (P.T.) the 
hexagonal (octagon) plaquette maps to the square (pentagon) 
plaquette. A local unitary transformation U is needed to 
bring the plaquette operators into the canonical form. 


The Hilbert space at site n is spanned by the occupation 
of the two fermion modes: |00)n, |01)n, |10)n, and |ll)n- 
The fermion parity of site n is given by 

Dn = (2V’L^an “ 1) (2V’^„V’/3n “ 1) = Inlnlnlt (21) 

The surface code model of spins can be obtained by pro¬ 
jecting to the subspace of even parity: = 1. This 

reduces the dimension of each site from four to two. 
We can map the allowed states {|00)n, |ll)n} to effec¬ 
tive spin states: = |00)n and |J|)n = |ll)n- Under 

this projection, the Hamiltonian described by Eqs. ( p^ - 
([T^ can be rewritten in terms of the effective spin 
Pauli operators. For instance, i7i7i|00)i = |ll)i/4 and 
i 7 « 7 ^|ll)i = |00)i/4, therefore ijiJi maps to S'f/4. It is 
straightforward to carry out this mapping and show that 
the surface code model of Majorana fermions indeed re¬ 
duces to the spin model in Eqs. @-(§ after projecting 
to the even-parity subspace. 

Eor the two unpaired Majorana modes shown in Eig.|^ 
the total fermion parity is given by P = '<^ 7 i 27 i 9 - How¬ 
ever, after the projection, P is no longer a physical opera¬ 
tor because it does not commute with the fermion parity 
operator P 12 and P 19 , namely it brings the system out 
of the even-parity subspace. A simple cure to this prob¬ 
lem is attaching a string operator 5 'i 2 ,i 9 connecting two 
unpaired Majorana modes. One possible choice of such 
a string operator is the product of all the link operators 
going from site 12 to site 19 

— ^7i27i9'^12,195 ^12,19 — ^12,11^11,10^10,9^9,20^20,19- 

( 22 ) 


This newly defined parity operator communtes with the 
fermion parity operator Using the mapping from 

Majorana operators to Pauli operators under even-parity 
constraint, Eq. (22) reduces to the expression of Eq. (16) 
obtained using Jordan-Wigner transformation. 


C. Perturbative limit of the Kitaev honeycomb 
model 

The surface code model can be obtained from the per¬ 
turbative limit of the Kitaev honeycomb model. It has 
been shown that the Kitaev honeycomb model admits 
unpaired Maj orana modes after introducing dislocations 
into the lattic^^^^. Naturally, one would expect the un¬ 
paired Majorana zero modes associated with twist defects 
in the toric code model can be deduced from the Kitaev 
honeycomb model with a certain type of dislocations. In¬ 
deed, we show this is the case by starting from the hon¬ 
eycomb model with dislocations and then performing the 
perturbation theory. 

The Kitaev honeycomb model h as sp in- 1/2 particles 
interacting through two-body terms^^EU 

Hhcm = -Jx a^<T^-Jy Y 

a^links ylinks zlinks 

(23) 

where the classification of x^y^z links is shown in 
Fig.||;a). This model can be exactly solved by mapping 
spin to four Majorana modes 7 ^’^’^ and 77 ^ 

< = (24) 

with the Majorana modes satisfying the anti¬ 

commutation relation 

{7m 5 7n} ~ 7n} — {7m’ 7n} = 0- 

(25) 

The dimension of the Hilbert space at each site is in¬ 
creased from two to four after the mapping. Such a prob¬ 
lem is fixed by defining a projection operator 

Dn = InllluVn, (26) 

and constraining the physical states to be eigenstates of 
Dn with eigenvalue +1. The Hamiltonian can be rewrit¬ 
ten in terms of Majorana modes 

^HCM = ^ X] X] (27) 

a=x,y,z (mn) 

where {mn) denotes the bond connecting spins m and n 
and the link operator H is straightforward 

to check that commutes with the Hamiltonian and 
hence can be set to either +1 or —1. The choice of the 
value of the link operator fixes the Z 2 background gauge. 
As a result, the Hamiltonian becomes quadratic in t erms 
of T] Majorana modes and can be solved exactl}!^^. 

A pair of twist dislocations can be added by modi¬ 
fying the bonds according to Eig.|^a). Every 7 Majo¬ 
rana mode is paired with another 7 in the link operator 
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except 7^2 and 7^9, both of which disappear from the 
Hamiltonian and hence are unpaired Majorana modes. 
Together they form a fermion mode introducing an addi¬ 
tional topological degeneracy of two. Next, we show that 
in the strong 2;—bond limit the honeycomb model with 
dislocations maps to the surface code model with twist 
defects, and the unpaired Majorana modes carry over. 

In the limit ^ > 0, it becomes energetically 

favorable to have the spins in the z—bonds aligned, i.e., 
— 1- Each 2;—bond becomes an effective spin with 
states 1^) and |J|) In this case, the 

Hamiltonian can be split into two parts 

^HCM = i^o + E, (28) 

Hq = —Jz (29) 

2;links 

V = -J,Y. « -JyH «■ ( 30 ) 

a^links ylinks 

Hq restricts the ground state to the subspace of |^) and 
|J|) and V can be treated as a perturbation. For a hexag¬ 
onal plaquette without dislocations as shown in Fig.j^b), 
the first non-vanishing term from the T^-order perturba¬ 
tion theorjEl is 

= T^VGoVGoVGoVT = - ^ + const., 

2; p 

( 31 ) 

where Go = {Eq — is the non-interacting Green 

function with Eq being the ground state energy of Hq, 
and T is the projection operator of the subspace |J| 
)}. Regrouping the Pauli operators in enables us to 
rewrite in the effective spin representation 


= -(a4Vf)KV|)(cT^2")(^f^!) ^ sys^.sysi 

(32) 

where 5'^’^ is the effective Pauli operator for the strong 
z-bonds shown in Fig.j^b). For a plaquette with a twist 
dislocation as shown in Fig.j^b), we follow the same pro¬ 
cedure and the 5^^-order perturbation theory gives rise 
to the following Hamiltonian in the effective spin repre¬ 
sentation 

2; p 

= sisisysisi. (33) 


A local unitary transformation U = will 

bring and into the same form as the surface 
code model introduced in Eqs. (i-©. Therefore, after 
the perturbation theory and the transformation U, the 
honeycomb model in Fig.j^a) reduces to the surface code 
model in Fig.[^ The unpaired Majorana modes 7^2 and 
7^9 associated with twist dislocations in the honeycomb 
model carry over to the surface code model. Similar to 
the case of Majorana fermion model in Sec. H B . ^7 i27i 9 is 



FIG. 4. (color online) Measurement of the parity (topological 
charge) of two twist defects in the surface code stabilizer set¬ 
ting. The direct way is to measure simultaneously the Pauli 
operators (red circles) making up the parity operator. An¬ 
other indirect approach is to create a measurement qubit of 
two holes by stopping two stabilizers shaded in dark brown. 
Moving one of the hole around the twists encodes the parity 
in the logical Zl oi the measurement qubit (product of the 
X’s connecting the two holes). By measuring the logical Zl 
before and after moving the hole, the parity of twists can be 
read out. 


not the physical parity operator because it does not com¬ 
mute with the projection operator D 12 and Dig. Again, 
the solution is to attach the link operator connect¬ 
ing site 12 with site 19 to ^712719 so that it commutes 
with Dn. In the perturbative limit, this newly defined 
parity operator reduces to the parity operator we found 
previously in surface code mode. It is highlighted in red 
circles in Fig.[2 


III. NON-ABELIAN STATISTICS OF 
MAJORANA FERMIONS 


In this section, we focus on demonstrating the non- 
Abelian statistics of Majorana fermions using the twist 
defects through braiding operations. Since introduc¬ 
ing the twists requires surgery of the underlying lat¬ 
tice Hamiltonian, moving twists would constantly change 
the Hamiltonian, which is not favorable from the 
point of view of experiments. Instea d, we propose to 
use measurement-based braidin^^^^^ to show the non- 
Abelian statistics. Our proposal requires only the par¬ 
ity measurements of a pair of twists. We first demon¬ 
strate how to perform such parity measurements in a 
surface code setting. Then we elaborate on the details 
of measurement-based braiding, and show that a single 
cycle of parity measurements followed by a single logi¬ 
cal qubit operation is sufficient to simulate the braiding 
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operation. This removes the uncertainty of number of 
measurements in the forced measurement schem^^. 


A. Parity measurement of twist defects 

All of our previous discussion is based on the Hamil¬ 
tonian formulation, and it is still very challenging to en¬ 
gineer the four-body or five-body plaquette operators in 
Eqs.([2])-@. Alternatively, one can use the stabilizer for¬ 
mulation in topological code^^^^^. In this case, there is 
no Hamiltonian and each term of the Hamiltonian is re¬ 
placed by a stabilizer operator. The stabilizers can be 
implemented using a combination of single-qubit gates 
and two-qubit CNOT gates followed by a projective mea¬ 
surement; see Fowler et al^ for details. In each cycle 
of the code, all the stabilizers are measured simultane¬ 
ously. The resulting state after a projective measurement 
is an eigenstate of all stabilizers with randomly selected 
eigenvalue of either +1 or —1. Error detection is done 
through repeating the cycle of stabilizer measurements, 
comparing the measurement results, and identifying spe¬ 
cific errors on particular qubits with the help of classical 
matching algorithms. 

To demonstrate the non-Abelian statistics of twists, we 
need to perform parity measurement frequently. We pro¬ 
pose to do such parity measurement in a stabilizer setting 
as shown in Fig.|^ The first method would be directly 
measuring the parity operator Otwists that is a string of 
Pauli operators highlighted in red circles in Fig.[^ This 
can be done in three steps. First, we turn off the sta¬ 
bilizers involving those Pauli operators to isolate them. 
Second, we perform corresponding Pauli XjYjZ mea¬ 
surement on individual qubits. The product of the mea¬ 
surement outcomes is the parity. Finally, we turn the 
stabilizers back on. Error detection is done by compar¬ 
ing the stabilizer measurements before and after turning 
off the stabilizers. 

A second indirect approach to measure the parity is to 
create an ancillary measurement logical qubit by stop¬ 
ping measuring two stabilizers shaded in dark brown in 
Fig.|4| This creates two holes, and generates a four¬ 
dimensional degree of freedom with each hole taking on 
an eigenvalue of +1 or —1. We effectively encode a sin¬ 
gle qubit in this four-dimensional subspace by defining 
the logical Zp as the operator connecting the two holes 
and logical Xl as the stabilizer of one of the two holes. 
Then, we move one of the holes (the top right one in 
Fig.|^ around the two twists forming a loop. This pro¬ 
cess performs a CNOT gate betwee n the ancillary qubit 
and the qubit encoded in two twist^^^^ 

Zl Zl 0 Otwists- (34) 

Therefore, by comparing the outcome of Z^ measurement 
before and after moving the hole around the twists, we 
can read out the qubit encoded in the two twists. Note 
that such a measurement is quantum non-demolitional 



FIG. 5. (color online) Three different bases of four Majorana 
modes, (a)-(c) Three different ways to group Majoranas 1 to 
4. (d)-(f) The corresponding fermion number states expressed 
in fusion tree diagrams. 


and can be repeated several times to improve the accu¬ 
racy of the parity readout. 


B. Basis transformation of Majorana fermions 


Now, we work out the details of measurement-based 
braiding by following the wave-function after each pro¬ 
jective measurement in the Schrodinger picture. Because 
the scheme involves measuring the topological charge of 
anyons in different bases, a basis transformation between 
different ways of grouping Majorana modes is needed. 
Given a set of four Majorana zero modes, we have three 
distinct ways to group them into two pairs with fermion 
charge i and j as shown in Fig.[^a)-(c). These ba¬ 
sis states can be represented using fusion tree diagrams 
shown in Fig.|^d)-(f). The Majorana fermions can be 
treated as Ising anyons (up to an overall irrelevant phase 
factor in the braiding operator). The Ising anyon model 
has three types of topological charges: vacuum /, Ising 
anyon a and fermion The fusion rule of Ising anyon 
modeP^ is given by 


I X I = I X a = I X 'ip = 

(J X (j = I -\r 'Ip^ crX'0 = cr, %p X %p = I. (35) 


Depending on the total fermion parity of the four Ma¬ 
joranas, the charges i and j further fuse into anyon 5 
which is either vacuum I for even parity or a fermion 
pj for odd parit3^^. The basis transformation can be 
viewed as changing the ordering of fusion. This amounts 
to transform the fusion tree diagrams in Fig.[^d)-(f) into 
each other. Such a transformation is dictated by the F- 












matrices, i?-matrices and 5-matrices as following 





where the F, R and B {= F ^RF) matrices for Ising 
anyons are given bjEH 



[Raa] = 


1 0 
0 i 






Here, the matrices are all in the basis of {/, t/^}. 

= [R^axi)]'y<^ = -1- All the other matrix ele- 
ments are either 1 if it is allowed by the fusion rules or 0 
if not allowed. 

We first work with the even-parity subspace which 
means anyon 5 after all the fusions is vacuum I in 
Fig.|^d)-(f). As an example, the steps to transform from 
the basis (13,24) to (12,34) are given by 


1 234 1 234 



Mathematically, this transformation can be expressed as 

\'i) 24: = [-^ 224 ] 2 7 [^ 123 ]*^ [-^/c34]i/c 1^)12 1^)345 (38) 

jk 

where \k)ab denotes the wave-function of the state after 
fusing anyons a and b into a channel with topological 
charge k. Here, 1—4 are Ising anyons a, j is a according 
to the fusion rules, and k can be either I or 'ip. Plugging 
in the matrices from Eq. (37) gives the following result 


|0)i3|0)24 

11)1311)24 


= 




(e) _ 

13^12 — 


e 8 

71 


13^12 

1 —i 
—i 1 


|0)l2|0)34 

11)1211)34 


FIG. 6. (color online) Measurement-based braiding involv¬ 
ing four Majorana modes, (a) The initialization of Majorana 
1 and 2 in the vacuum state, (b)-(d) Measurements of the 
fermion number of Majorana 1 and 3, 1 and 4, and 1 and 2, 
respectively, (e) Definition of the single-qubit logical A, Y 
and Z operators. 


Repeat the above procedure, we obtain the unitary trans¬ 
formations between (14,23) and (12,34), and between 
(13,24) and (14,23) in the even-parity subspace 


^14^12 - ^ 


1 1 

—i i 


U- 


(e) _ 

13^14 — 


e-*8 h -I 


(40) 


Similarly, the basis transformation in the odd-parity 
subspace is obtained in the same way as 


|0)i3|1)24 

11)1310)24 


— ^13^12 


^13^12 


^14^12 


e 8 


5, 


(o) 

13^14 


V 2 
1 

71 


1 -i 

1 


72 


i —i 

1 1 

1 

-1 


|0)i2|1)34 

1210)34 


(41) 


It is straightforward to check that the unitary transfor¬ 
mations satisfy the consistency equation 


^ 13^12 — ^ 13 ^ 14 ^ 14^125 


Tj{o) _ tjKO) jjy 

^ 13^12 ~ ^ 13 ^ 14 ^ 14 ^ 12 - 


Ao) 


jio) 


( 39 ) 


( 42 ) 
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C. Measurement-based braiding 


pending on the outcome the state is 


With the unitary transformations between bases in 
both the even-parity and odd-parity subspaces, we 
are ready to work out the details of measurement- 
based braiding (MBB) in the Schrodinger picture. The 
measure-only approach to topological quantum compu¬ 
tation was proposed by Bonderson et alJ^ to imple¬ 
ment the topological gates without physically braiding 
the computational anyons. The braiding operation in¬ 
stead is replaced by a series of quantum non-demolitional 
topological charge measurements as shown in Fig.^a)- 
(d). 

We consider the special case of Ising anyons (Majo- 
rana fermions) only. In our case, the topological charge 
(fermion-parity) of each pair of Majorana modes can 
be measured using the scheme proposed in Sec. [mg 
Anyons 3 and 4 are the computational anyons to be ex¬ 
changed, and 1 and 2 are auxiliary anyons employed to 
assist MBB. The scheme works in the following way: one 
first initialize the two anyons 1 and 2 in the vacuum state 
with n ^2 = 0 as shown in Fig.j^a). Then a forced mea¬ 
surement is perform on anyons 1 and 3: if riis is not 0, 
we go back to measure ni 2 followed by nis and repeat 
until we obtain nis = 0 [Fig.j^b)]. In some sense, we 
force anyons 1 and 3 to fuse into the vacuum sector by 
repeated measurements. Physically, this procedure real¬ 
izes anyonic teleportation of the state encoded in anyon 
3 to anyon 2. Similarly, forced measurements are done 
on anyons 1 and 4, and 1 and 2 shown in Fig.[^c)-(d) 

to have ni 4 = 0 and = 0? respectively. It turns out 
the resulting effect is equivalent to a braiding operation 
on anyons 3 and 4^. Due to the inherent probabilis¬ 
tic nature of forced measurements, the operation time 
of each measurement-based braiding unavoidably varies 
from run to run. This imposes difficulties of synchroniz¬ 
ing the clock if one wants to use MBB to perform topo¬ 
logical gates. Recently, the generic case of all possible 
intermediate measurement results is considered^. Here, 
it is our aim to study the generic case in detail and we 
find that one can perform MBB with a fixed number of 
measurements removing the uncertainty associated with 
forced measurements. 

We assume that after the initialization step in Fig.j^a), 
the state of the four anyons is described by 

|^(^)) = |0)i2 [(a|0)34 + ^11)34] • (43) 

Using the unitary transformations between bases in 
Eqs. ([3^-([4T]), we can rewritten the initial state in the 
basis of (13, 24) 

p—i7r/8 

= ^[«(|0)l3|0)24 + i|l)l3|l)24) 

+/^(|0)i3|1)24 + '<^ 11 ) 1310 ) 24 )] • (44) 
After the projective measurement of ni 3 in Fig.j^b), de- 


|^(^)) 


e "8 |0)i3[a|0)24+ ^ 11 ) 24 ], ifni3 = 0 
[q^| 1)24 + /5|0)24] , if ni3 = 1 


(45) 

Next, we rewrite |T^^^) in the basis of (14,23) and per¬ 
form projective measurement of ni 4 shown in Fig.j^c). 
The resulting state is 


{ e *4 |0)i4[a|0)23+^ 11 ) 23 ], ni3 = ni4 = 0 
|l)i4[ - a|l)23 + ^ 10 ) 23 ],^13 = ni4 = 0 
ie“* 4 |0)i4 [(a|0)23 — /^|1)23], ^13 = f^i4 = 1 

|l)i4 [q^|1)23 + ^ 10 ) 23 ], rii3 = ni4 = 1 

(46) 

where n = 1 — n. Finally, we transform into the 

( f 

basis of (12,34) and carry out the measurement of n \2 
in Fig.j^d). Up to an overall phase factor, the final state 
can be grouped as 




|0)i2 [0^10)34 + i/3|l)34], ni3 = ni4, = 0 

|0)l2 [q;|0)34 - i/3|l)34] , ^13 ¥= '”-14, ^12^ = 0 

|l)l2 [*Q;|1)34 + /3|0)34] ) ^13 = ni4, = 1 

Jl)i2 [iQ!|l)34 - /3|0)34] , ni3 ^ nu, = 1 


(47) 

A close look at the final state tells us that if we apply an 
operator P as defined below to the final state, it is the 
result of a braiding operation on the initial state 




(48) 


where Rza = (1 + 7473)/\/2 and 73,4 are the Majorana 
operators of anyons 3 and 4. The operator P depends on 
the measurement outcomes and is given by 


[, 

if) 

ni3 = ni4,nF2 

= 0 

Z = i7374, 

ni3 ^ ni4,n[{'' 

= 0 

Y = ^7174, 

If) 

ni3 = ni4,nf2 

= 1 

X = i'jijz, 

^13 7^ ni4,n^'^^ 

= 1 


(49) 


Operator P has an intuitive interpretation: if we encode 
information in the parity of anyons 3 and 4, then P is the 
single-qubit operator manipulating the degree of freedom 
associated with anyons 3 and 4 as shown in Fig.j^e). 
By monitoring the measurement outcomes, we apply one 
of the logical Pauli operators from the set {/, X,T,Z} 
to complete the measurement-based braiding in Eq. ( [4^ . 
In the surface code setting, X, Y and Z are the logical 
operators connecting the twist defects which commute 
with the stabilizers. Eor example, Z operator is defined 
as a string of Pauli operators highlighted in red circles in 
Eig.[^and Eig.|^ 

We can envision using the MBB to demonstrate the 
non-Abelian statistics of Majorana fermions. A minimal 
set of 6 twist defects would be enough. In addition to 
twist 1 to 4 in Eig.[^ we have to include twists 5 and 6. 
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We can initialize in the basis of (35,46) with the vacuum 
state. Then, we perform a MBB to braid 3 and 4. If we 
measure the parity of 3 and 5: P 35 , there will be 50% 
probability that the parity is changed. Generally, we 
can perform n MBB of 3 and 4, and then measure P 35 . 
The probability of parity change will be 1, 0 for 

n = 1, 2,3, 0 (mod 4]P^. This is an appealing alternative 
to the condensed matter settings because the necessary 
experimental capacities are within reach due to recent 
advancements in the surface cod^^^M^. 


IV. CONCLUSIONS 

In conclusion, we have revisited the surface code model 
with twist defects. Using a variety of different techniques, 
we have shown explicitly the emergence of unpaired Ma- 
jorana zero modes associated with twist defects. We also 
propose a scheme to measurement the parity (topological 
charge) of pairs of Major anas. The parity measurement 


serves as a building block for measurement-based braid¬ 
ing. We investigate the possibility of performing such 
braiding without forced measurements. It turns out that 
it can be done with a cycle of three topological charge 
measurements and one additional logical operation. The 
uncertainty associated with forced measurements is re¬ 
moved by our approach. This makes measurement-based 
braiding an appealing method for both demonstrating 
non-Abelian statistics of Majorana fermions and build¬ 
ing Clifford gates for topological quantum computation. 
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